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introduction 


The public sector affects both the level and 
distribution of national income. Until fairly recently 
concern has been focused largely on questions of the 
level. This has been reflected in Canada in the development 
of a number of aggregate economic models which have 
tended to concentrate on the business sector. Perhaps 
the best known of these are RDX2 and CANDIDE. Although 
these models do provide treatment of both the household 
and government sectors, the manner in which this is 
done does not lend them to the study of the income 
distributional consequences of government programs. 
This is the case because the distribution of income is 
most meaningfully considered in relation to individual 
households or families, the recipients of national 
income, and these models do not maintain sufficient 
household or family detail. It is clear that any 
exercise designed to elucidate distributional issues 
must begin with the household sector portrayed in great 


particularity, preferably at the level of the individual 


person. 


Of course, given the availability of elaborate 
sets of microdata, it is possible to develop fairly 
simple models that enable one to pose "what if" questions 
for some past period. But this kind of static analysis, 
although useful, does not allow one to come to grips 
with a host of questions bearing on income distributional 


issues which are crucially dependent on time. To be 
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relevant, it is also necessary to have the ability to 
forecast the workings of a large number of factors, 
economic and demographic, that bear on the household 


sector and determine the distribution of income. 


It is also necessary to have the capacity to model 
government programs as comprehensively as possible, and 
in considerable structural detail. It is not possible 
to know what any single program will do if we are 
unable to situate that program in the context of an 
environment produced, in part, by a number of other 
programs. Further, it is not enough to model the 
general direction of the effects produced by a given 
program as a whole. We are concerned to understand the 
Significance of particular program designs. To do this 
we require the ability to test the effect caused by an 


alteration of internal program components. 


The term simulation is usually used to describe 
techniques related to the construction of models or 
simulators whose operations are intended to resemble 
the behaviour of actual or potential operating systems. 
Microdata simulation involves the construction of 
simulators intended to function in the same manner as 
operating systems comprised of a large number of basic 
components or decision units. In carrying out this 
form of simulation one may employ the technique of 
endowing microcomponents with behavioural functions and 


of deriving the consequences (in terms of individual 
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behaviour) of different environments (specified by a 
set of independent variables) on the microcomponents. 
Or, One may employ the technique of stochastic events. 
In this case microcomponents are not endowed with 
explicit specifications as to their behaviour under 
differing environments but rather are thought to behave 
as particles governed by specified laws of chance. 

That is, the causal relationships that actually determine 
behaviour are considered only implicitly, either 
because they are too complex to be modelled, or because 
insufficient data exists to estimate the specified 
model. One is only able to observe that the micro- 
components reside in a certain "state" for a period of 
time and then move to other "States". The precise laws 
governing individual movement are unknown. What is 
known is that if a large number of "State" changes or 
movements are observed, the behaviour in question can 


be described as if it depended solely on chance. 


The probabilistic analog of a deterministic causal 
process is the first order Markov-chain. In this kind 
of process, the future state of the world is completely 
independent of time periods preceding the present. In 
the present context, then, micro-component behaviour 
would be described completely by its present state, and 
the Markov-chain probabilities that relate the present 
to the future. In the broadest possible terms this is 
the POLSIM Model: it is a Markov-chain simulation of 
individual demographic, labour force, and market income 


behaviour. 
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More specifically, POLSIM is an annual microdata 
model of the Canadian household sector. The basic 
component of the model is the individual Penconn 
Individual persons may be associated into nuclear 
family units in the model but the prime focus is always 
Maintained on the individual. The model receives, as 
exogenous input, a specification of the native Canadian 
population for some year. This specification is made 
in terms of a number of characteristics (a state vector) 


for every individual in the initial population. 


The individual state vector is comprised of three 
different sets of characteristics: demographic (e.g., 
age, sex, etc.), activity (e.g., weeks employed, weeks 
unemployed, etc.) and income (e.g., annual wages, 
annual dividend income, etc.). Most of these character- 
istics change over time. The function of the model is 
to effect these changes in the light of the individual's 
particular circumstances, as portrayed by his state 
vector immediately before the change, and conditions 


prevailing in the socio-economic environment. 


The POLSIM model is constructed as a number of 
connected blocks: (i) immigration, (ii) demographic 
(iii) activity status, and (iv) market income. Each 
of these blocks contains a series of processes or 
transformations which operate on individual state 


vectors to produce annual change. The structure of the 
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model is given in Figure 1. The model begins with a 
description of the population at some year t, and then 
proceeds to synthesize a new population representative 
of the year ttl. This updated population can then 
become input to a second run of the model, and so on, 
until the initial population has been updated as many 


years as required. This final population (as well as 


any intermediate population) can then be used to simulate 


the absolute and distributional consequences of a given 


government program. 


The next section of this paper will describe the 
individual state vector and the base year population in 
more detail, while the third section will elaborate on 
the model structure. Section 4 will then outline some 


simulation results. 


The POLSIM Data Base 


The data base for the POLSIM model is obtained 
from the Statistics Canada Survey of Consumer Finance. 
In the recent past this Survey has been conducted every 
two years. In future, it will be taken annually, but 
every second survey will be of small scale and of a 
specialized nature so as to make it unsuitable for 
initial year model population specification. For the 
present purpose, therefore, we may only regard each 


second survey as useful. 
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The SCF is at present the most extensive micro-data 
base available in Canada, with the exception of the 
census. Coverage extends to all of the population of 
Canada with the exception of (i) persons resident in 
the territories, (ii) persons resident on Indian reservations, 
and (iii) persons resident in institutions (e.g., 
prisons, mental hospitals, etc.). The Survey carries a 
wide variety of income and other information capable of 
organization on either an economic or census family 
basis. The quality of this data is generally very high 
and its availability (for our present purpose) every 


two years makes it timely. 


In the last section we described in very broad 
terms the nature of the individual state vector, that 
is, the fact that it contains demographic, activity and 
income characteristics. The complete vector in its 
most detailed form contains 23 characteristics in all. 
These characteristics are derived from the individual 


SCF records and are listed below. 


Ales Previous Year Family Unit 

2 Family Unit Identifier 

3. Province 

4. Size of Family 

Dn Census Family Relationship 

6. Marital Status 

gs Age 

8. Sex 

o. Major Source of Income 

10. Weeks in School 
Le Weeks employed 

iZs Weeks unemployed 
ESBS Weeks in the non-labour force 
14. Education 

LD April Activity Status 

16. Weight 
kis Employment Category 
athe Employment Income 
zo, Interest and Other Investment Income 
20; Dividends 

2h5 Retirement Pension, Superannuation, and 

Annuities 
eae ae Other Money Income 


255 Total Income 
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These 23 characteristics taken together "describe" 

one individual. The entire Canadian population could 

in principle be "described" in the context of the model 

by a large number of individual state vectors - one for 
each individual in the population. This would, however, 
be a very inefficient procedure. A much more economical 
means of accounting for the entire population is by way 

of a representative sample. In this case, each individual 
state vector in the model stands for or represents a 


larger number of "identical" persons in the real population. 


The model requires, as initial input, a weighted 
sample of individual state vectors which describe the 
population of Canada for some year. At the present 
time we have been working with two such base years, 
1967 and 1971. The Survey of Consumer Finance undertakes 
its geographically stratified sampling of Canadian 
households in April of the year following the survey 
year (e.g., April 1968: tor 1967 or April 1972 for 
1971). The sample included 37,985 individuals over 
fourteen years of age and in receipt of cash income in 
1967 and 43,039 individuals aged fourteen years or more 
who were in receipt of cash income in 1971. Over time, 
of course, aS more current data becomes available, 
initial year state descriptions for more recent periods 
can be constructed. It is important to be able to 
continually update the initial year as a guard against 


the generation of inaccurate projections. 
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Individual state vectors are weighted by the 
Survey of Consumer Finance in order to calculate population 
estimates. These weights are unequal, running from 30 
up to 3,000 in increments of 10. In order to conform 
to the logic of POLSIM, individual state vectors must 
possess equal weights. It was therefore necessary to 
adjust the SCF samples in such a way as to produce 
equally weighted records. We adopted a common weight 
of 50 (i.e. one individual in the sample represents 50 
in the real population). The adjustment was then done 
by an unbiased rounding procedure (rounding to multiples 


of 50) combined with record replication. 


The process of replicating sample state vectors 
increases the number of sample individuals from 79,479 


to 397, 960" in 1967 ‘and “from /9.,528° to 425,864 an- 1971. 


The Model Structure 


As mentioned previously, the POLSIM model is 
divided into 4 blocks: Immigration, Demographic, 
Activity, and Market Income. The purpose of the 
Immigration Block is to fabricate the state vectors of 
the immigrant population that will arrive in Canada 
during the year being simulated. The purpose of the 
other three blocks is to update respective subsets of 
the individual state vector. This section will briefly 
describe these four blocks, outlining their functions 
and the way in which they operate. The reader is 
referred to the Appendix for the detailed flow chart of 


each block. 
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The Immigration Block 


The objective of the immigration block is to increment 
the base population in a given year by the total number 
of immigrants that will arrive in that year. This 
entails two distinct problems: the projection of the 
total number of immigrants who are to arrive in the 
given year, and the synthesis of individual state 
vectors for each of these individuals. The present 
version of the model does not attempt to determine 
projections of the total number of immigrants. Such 
projections are simply assumed by the model to be 
exogenously determined. The immigration block is 
therefore concerned principally with the construction 
of individual state vectors for a pre-determined number 


of people. 


The immigration block operates as follows. The first 
step is to calculate the total number of people in each 
of 240 age-sex-marital status-province classes. These 
classes, consist only of adults. The numbers of children 
in given classes are determined after married couples 

are formed. The number of adults in a given class is 

the product of the total number of immigrants arriving 
and the probability of being in a given class, divided 


by the weighting factor of 50. 
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The model then calculates the individual state vectors 
on a province by province basis. State vectors are 
determined for all the people in a given province, and 
these are then written out on tape. Those for single 
males are calculated first. The procedure is to consider 
one individual at a time, calculating his entire state 
vector from various probability distributions. The 
next individual is then considered, and so on, until 
all single males have been dealt with. An analogous 
procedure is then repeated for all single females, and 
finally for all members of families. The family loop 
includes the calculation of the number of children a 
given family will have, as well as these children's 


respective state vectors. 


Once all the immigrants in a given province have been 
created, their state vectors are written out, and the 
calculations are repeated for the next province. The 
final output from the immigration block consists of 
state vectors describing all of the immigrants in each 
of the ten provinces. These state vectors are now 
ready to enter the demographic block along with the 


remainder (i.e. non-immigrant portion) of the population. 


The Demographic Block 


The purpose of the Demographic block is to update the 
demographic variables of the individual state vector. 


These variables are the individual's province of residence, 
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his family size, his dependency status (whether he 

is a family head, a wife, or a dependent child), his 
marital status, and his age. The Demographic Block 
also determines whether or not a family will emigrate. 
Since the basic time unit of the POLSIM model is one 
year, the updating of these state variables is on an 


annual basis. 


There are three features to the Demographic Block that 
distinguish it from the other three blocks in POLSIM. 

The first is that it eliminates and creates individual 
records. New records are created in the Demographic 
Block through the simulation of births, while others 

are eliminated through the processes of death and 
emigration. The second distinctive thing about the 
Demographic Block is that it proceeds in two distinct 
phases. The Activity and Market Income blocks are 
somewhat simpler in this respect. They require only 

a single sequential processing of all individual records. 
Two phases are necessary in the Demographic Block 
because of the need to simulate marriage. Since 
individuals in the base year population must marry 

other individuals in the same population, it is 
necessary to first determine all of the "marriageable" 
individuals. Only then can these individuals be paired. 
The final distinctive feature of the Demographic Block 

is that if deals with family records, whereas the other 
two processing blocks work with individuals only. The 
reason for this is that demographic variables are inherently 
familiy variables; that is, dependency status, family size 
etc., can clearly only be determined within the context 


of a whole family. 
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The Demographic Block consists of two distinct phases. 

The input to the first phase is the file containing 

the native Canadian population of the previous year 
together with the population of new immigrants for 

the current year. The first phase of the block then 
produces three output files. The first output file 
contains the records of all individuals who either 
emigrate or die during the year being simulated (i.e., 

the emigration and death registries). The second output 
file contains the records of all individuals who have 

been determined to be eligible for marriage during the 
simulation year. This file, called the "Marriage Pool", 
will be used in the second phase of the Demographic 

Block. The third output file contains the output records 
of all remaining individuals (all those who have not 

died, emigrated, or been deemed ready for marriage). 

In phase one of the Demographic Block, then, all individuals 
who have not died or emigrated will have had all their 
demographic variables completely updated. Persons who 
will marry, on the other hand, are updated in all respects 


except for "pairing" and their province of residence. 


The second phase is exclusively concerned with persons 

in the marriage pool and determines "who will marry 
whom". The individuals who have previously been recorded 
in the marriage pool (i.e., the second file mentioned 
above) are formed into couples on the basis of age, 


province of residence and education. The demographic 
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characteristics of these couples, with the exception 

of province of residence, have all been updated in 
phase one, and phase two now determines the province 

of residence of the new couples, in the same manner 

as is true of persons who are not marrying between 

t and t+l. All the couples in this "Marriage Pool" 
file are then merged with the records in the third 
output file (i.e., individuals who do not die, emigrate 


or marry during the simulation period) of phase one. 


The basic method by which the stochastic decisions 
pertinent to the updating of the state vectors are made 
is Monte Carlo simulation. The first decision to be 
made after a family record has been read is that of 
whether or not the individuals in the family will survive. 
Survival is assumed to be an individual matter. All 
individuals within the given family are processed one 
by one. If the Monte Carlo procedure determines that 
the person shall die, he is recorded on the death 
registry, and his record is not considered further. If 
nobody in the family survives, another family record is 
read from the Initial Year file. Surviving individuals 


proceed to the birth process. 


Births may be legitimate or illegitimate. If the woman 
is married, and it is determined that she will have 

a child, the child is simply added to the given family 
unit. If, however, she is not married and she is not 
the head of a family unit, she and her baby are assumed 
to form a new family unit; and, of course, she becomes 


the head of this new unit. 
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The divorce process is now commenced for all families 
containing both a head and a spouse. If divorce is 

determined ‘to’ occur ,) the family as*split in ‘two. All 
children are assumed to go with the mother into a new, 


separate family unit with the mother as head. 


All non-married individuals, excluding those who have 
become divorced in the current simulation, are now 
tested to determined whether they will get married. If 
this event occurs, the person is declared to be the 


head of a new family unit. 


All individuals in each of the original families who 
are dependent (i.e., neither heads nor wives) are now 
tested to see whether they will leave home or not. 

If this event occurs, the individual is declared to 

be the head of a new family consisting of himself 
alone. It is clear that each of the original family 
units might generate a whole set of new family units 
through marriage of dependents, divorce, birth of an 
illegitimate child, or as a consequence of a dependent 


leaving home. 


For all families it is necessary to determine whether 
or not they will emigrate. This is considered to be 
a family decision, rather than an individual decision. 
All families that do emigrate are recorded on the 


emigration registry, and are not considered further. 
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Those family units that do not emigrate, with the 
exception of those who are going to be married in the 
current simulation, are now processed to determine 
their new province of residence. As mentioned above, 
province of residence of newly married couples is 
determined in phase two, after marriage has taken 


place. 


The final step of the first phase is to record all 

family units on the proper files. Those families whose 
head is to be married are recorded on the marriage pool 
tape (file #2), while all others are recorded on file 

#3. (It will be recalled that file #1 contains the 

death and emigration registries). At this stage file 

#3 contains all remaining families, with the exception 

of new families formed through marriage. All demographic 
variables of the individuals on the third file have 


been completely updated. 


In the second phase of the Demographic Block the individuals 
who are to be married are processed to form couples. 

Once this is completed, all the new husbands and wives, 
together with any dependents they might have, are 

formed into family units. The province of residence of 
these families of file #3 is then stochastically determined 
to create the completely updated initial year file. 

This is the final product of the Demographic Block, 

which now becomes input to the Activity Block that 


follows. 
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The Activity Block 


Broadly speaking, the purpose of the Activity Status 

Block is twofold: first, to update those variables 

in the individual state vector that describe what 

a person is "doing" during the year being simulated; 

and, second, to make certain adjustments to the individual 
state vector preparatory to updating the person's income 
in the Market Income Block. Specifically, the Activity 


Block does the following: 


(a) it determines the number of weeks during the 
year being simulated that the person spends 
in school, employed, unemployed, and in the 
non-labour force (where "non-labour force" 


is defined so as to exclude those in school); 


(b) it determines whether a person advances his 


education; 


(c) it determines the person's activity at the end 
of the year being simulated or the beginning 
of the next simulation year (the "year" being 


simulated is defined as April through March) ; 


(ad) it makes any necessary changes in a person's 


employment category; 


(e) and finally, it converts a person's annual wages 
to a weekly wage rate if he was employed as a 
Class B person in the year preceeding the one 


being simulated. 
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(A “Class B" person is one who is subject to 
unemployment). The Activity Block makes the requisite 
changes in the individual's state vector on the 

basis of his present year's demographic characteristics, 
his past year's activities, and exogenously determined 


monthly Canadian aggregate rates of unemployment. 


The methodological approach taken by the Activity Block 
is somewhat different than that taken by the other 
blocks in POLSIM. In the Market Income Block, for 
example, the various components of an individual's 
income are simulated directly. A transition matrix 
determines to what extent, if any, a particular 
component of a person's income is to be increased or 
decreased during the simulated year. In the Activity 
Block such direct transitions are not possible. Instead, 
a person is thought of as being in one of four activity 
states: employment, unemployment, non-labour force, 

or school. Transitions among these four states then 
take place month by month, tracing out a year long 
"activity history" for a given individual. The 

number of weeks in each of the four states, any change 
in education status, and any changes in employment 


category are then inferred from this history. 


Much of the Activity Block is concerned with the problem 
of calculating the transition matrices that are used 
to determine how a person will move from month to month. 
That is, these matrices are not read in as data, but 


rather are calculated from other data in the model. 
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A general picture of the Activity Block can be obtained 
from the Macro Flow Chart given in the Appendix. The 
block begins by reading in all of the input parameters. 
These parameters include regression coefficients and 
unemployment rates from which an initial set of transition 
Matrices is calculated, parameters that are used to 

adjust these probabilities in various ways, and probabilities 
used to determine whether and how a person advances 
through school. Once all of the parameters are properly 
constructed, individuals are passed through the logic 

of the simulation model proper. In this section persons 
are directed through various processes, depending on 

how they relate to the labour force (their Employment 
Category*), and a few other factors. Children, for 
example, bypass most of the Activity processes completely. 
And persons who are assumed to never become unemployed 
bypass the monthly labour force simulation. The relevant 
state variables are then updated, and the whole procedure 


is repeated for the next individual. 


The idea behind the whole Activity Block procedure is 
to infer what an individual does during a year by 
examining what he does in every month of that year. 

A year is assumed to run from April through April (thus 


providing a starting point for the next year's simulation). 


The employment category indicates whether this is an 
individual's first year in employment; whether he is 
retired and if it is his first year of retirement; 
whether he is in the labour force or not; and if he 
is in the labour force whether he will receive a 
private pension on retirement, and whether he will 
be subject to unemployment. 
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An individual thus starts out in April doing something 
(being in a particular education state, being employed, 
etc.) and this then influences (together with the 
exogenously inputted influence of aggregate unemployment) 
what he will be doing in May, and so on. The whole 

year is traced out in this way. An individual moves 
from state to state, depending on his present state, 

his age, the month being considered, and other factors. 
Once the entire year has been completed, the year 


itself can be summarized in the relevant state variables. 


The Market Income Block 


The POLSIM model divides an individual's total annual 
income into four main components: employment income, 
property income, retirement income, and other money 
income. Employment income is income in the form of 
wages and salaries, military pay and allowances, and 
net income from farming, fishing, and other forms of 
self-employment. Property income is divided into two 
subcomponents: dividends and other investment income. 
Dividends are self-explanatory. Other investment 
income consists, generally, of income from fixed face 
value assets such as bonds, deposits, and savings 
certificates, and all other forms of investment income. 
Retirement income consists of private pensions, super- 
annuation, and annuities. Other money income includes 
income from roomers and boarders, alimony, gifts, and 


income from any other source not mentioned above. 
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The purpose of the Market Income Block is to update 

these income variables, on an annual basis, as the 

person moves from year to year through the simulation. 
The main technique for doing this is again Monte Carlo 
Simulator. Broadly speaking, the updating process 
consists of two general problems. The first is to 

assign initial component incomes to a person if he is 
"eligible" and if the particular component being examined 
was zero in the previous year. (A "component" of 

income is one of the four sources mentioned above.) 

The second problem is to determine transitions on each 

of the income components that were not zero in the 
previous year. Both of these two kinds of processes 

are generally handled by either a deterministic function, 
a stochastic function, or a combination of these. The 
exact meaning of "eligibility" in the various cases, 

and the way in which the several transitions are effected, 


will now be discussed. 


The Market Income block begins by reading and assembling 
all of the parameters that will be necessary to update 

an individual's income variables. All of the constant 
parameters, which consist of transition matrices, 

initial income arrays, growth factors, etc., are read 

in from a single data tape. They are then simply 

stored for use by the relevant subroutines. The variable 
parameters are the particular year being simulated and 
the rate of inflation that is assumed to apply throughout 
the simulated period. Both of these variables are read 


in from cards. 
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After reading in the rate of inflation, the program is 
able to calculate a consumer price index vector. This 
is simply the consumer price index for the 15 years 
from 1967 to 1981. It will be used to calculate money 
growth in certain components of a person's income. The 
CPI vector is the last parameter necessary to effect 
individual transitions. All of the other sets of data 
have already been stored, and the program is able to 


begin to read and process individual records. 


Before the income updating processes begin, the program 
checks to see if the individual is less than 14 years 

of age. Children less than 14 are assumed by POLSIM to 
receive no income, and consequently the income variables 


of all such children are not updated. 


All other persons are assumed to be eligible for employ- 
ment, property and other kinds of income. For employment 
and retirement incomes, which depend on one's relation 

to the labour force, a variety of processes exist. The 
exact transformations that these two sources of income 
are subject to depend on certain elements in the indivi- 
dual state vector. That is, they depend on just what 


"kind" of person the given individual is. 


The main criterion for determining how a person's 
employment and retirement income will be updated is his 
employment category, which, as mentioned earlier, 


defines the way in which he relates to the labour 
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force. If he is a person who has entered the labour 
force for the first time in the current year then he 
must have an initial income assigned to him. His 
retirement income will be zero. If the person is a 
student then an employment income must be assigned if 
the student worked during the summer. If the person is 
a non-retired member of the non-labour force then no 
employment income is assigned. In neither of the above 
two cases is a retirement income assigned. If a person 
is retired, then his retirement income is assumed to be 
the same as in the previous year, and his employment 
income is zero. If a person is in the labour force and 
is assumed to be employed for the full year, then his 
employment income is determined by an annual employment 
income transition. His retirement income is of course 
zero. If the person is in the labour force but is 
assumed to be subject to unemployment then the Activity 
block has previously determined how many weeks he has 
worked, and the Market Income block now determines any 
change in his weekly wage. His employment income for 
the year is then the product of these two variables. 
His retirement income is zero. If the person is in his 
first year of retirement (he has just reached age 65), 
and he is not eligible for a private pension, then in 
the present version of the model he is assumed to have 
no private retirement or employment income at all. If 
a person has just retired and is assumed to be eligible 
for a private pension, then his pension is calculated 
and becomes the whole of his retirement income. His 


employment income is zero. 
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It will be recalled that property income consists of 

two components, dividends and other investment income. 

The transition processes that apply to dividends are 
exactly equivalent to those that apply to other investment 
income. The only difference between the two is that 

the probabilities defining the two processes are different. 
But the processes themselves, the definition of the 

cells in the transition matrices, and the variables 

on which the matrices are stratified are exactly the 


same. 


The updating of "property income" can be divided into 

two sections, depending on how much initial property 
income the person in question has. The first applies 
only to persons whose property income is less than $250 
and embraces those persons who are moving from zero to 
non-zero property incomes for the first time. For 

people in this group, transitions to other classes 

depend on both age and income. A random number samples 

a cumulative probability distribution which is contingent 
on the person's age and income. The sampling defines 

the new property income class, and hence the new property 


income (which is taken to be the midpoint of the class). 


The second set of transitions, which apply only to 
those whose initial property income is greater than 
$250, is virtually identical to the first. The only 


difference is that the transition matrices depend on 
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age only, rather than age and income. The idea behind 
this distinction is that total income is relevant in 
determining the level of property income a person is 
likely to achieve in the first instance. That is, the 
larger the person's total income, the larger his savings 
are likely to be. And the larger the person's savings, 
the larger his initial property income. Whether this 
initial amount of savings is then further built up or 
drawn down depends mainly on the stage of the life 


cycle, i.e., age. 


Other money income consists of room and board income, 
alimony, and other small items of income which are 
difficult to simulate. Because this kind of income is 
small, and confined to very few people, it is handled 

in a simple deterministic way. The only persons allowed 
to have other money income are those who had some in 

the previous year. The size of the current year amount 
is increased, however, to reflect any changes in the 


CPI. 


Once all of the components of an individual's income 
have been determined, his total income for the simulated 
year is computed by simple summation. This total is 


then recorded in the individual's state vector. 
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Model Simulation: 19 6 to. TOF i 


The complete POLSIM model was tested by simulating 
the four year period from 1967 to 1971. The initial 
year population for the simulation consisted of 397,960 
individual records derived from the 1967 Survey of 
Consumer Finance. The input parameters to the simulation 
(unemployment rates, total number of immigrants, and 
rates of wage inflation) were the actual values that 
obtained over the four year period. The present section 


discusses the results of this simulation. 


The results of the simulation of new immigrants 
over the four year period are summarized in Table l. 
It can be seen that the model slightly underestimates 
the total number of immigrants, even though the actual 
number is an exogenous input. The reason for this is 
that children are created stochastically in the Immigration 
Block, and that slight adjustments have to be made to 
the number of married women so as to equate them with 


the number of married men. 


The results of the Demographic Block processes are 
summarized in Tables 2 to 5. The aggregate totals 
produced by the Demographic Block are compared with 
actual figures derived from Vital Statistics. In 
general, the model performs very well. The tables 
indicate that the model tends to underestimate the 
number of live births, deaths and divorces while overestimating 


the number of marriages. 
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The results of the birth process are given in table 

2, and it can be seen that the simulation begins with 

a fairly large error in 1968 which then declines until 
1971 when the results are almost perfect. Part of the 
error in the early years (approximately 15,000 live 
births) can be attributed to initial population errors. 
In particular, the initial population is largely underestimated 
for women in the 20-24 age group, and since this is the 
prime child-bearing age, a large underestimate in 

births will naturally ensue. The underestimate in 
births declines as the simulation progresses for two 
reasons. First the model uses stationary fertility 
probabilities estimated for the latest year. The birth 
rate declined over the period 1967-1971, and since this 
was not reflected in the probabilities, one would 

expect any underestimate to decline as time goes forward. 
Second, the underestimate will also decline as time 
progresses because the effect of the initial population 
error will decrease. In the initial population, females 
in the 15-19 age group are only slightly underestimated 
(as compared with the large underestimate in the 20-24 
age group). Therefore in each succeeding year, the 20- 
24 female age group will more closely approximate the 
true population in that age group, and hence the underestimate 
in births which results from an underestimate in the 
20-24 population will progressively be eliminated. In 
examining the regional results, table 2 indicates that 
there is no evidence that fertility probabilities 


should be regionalized. 


motuimnteeebhea! Yeuies wt auiae 


ond el etuz ebnte. ona: 
as laialibdiaiied 


an 
td pit dims veewastiate et euene ytteugdlien tle wth 


er 48 ‘ab eapapertg jotdutoiehe we’ ns asniivel erigxte 


erilisze} ywewoldess aden Intom past rast ea 


ajzid sft .sesy Je53nf ond Ybe Busamisee ealtilidsdesg 
eit aSaty oom ~LVCIWTOEL Potueg mae ewe Peewee 
hisov oi 2929 L0ideciong 923) i) Saguetien Som eee 


fagutc’ Heep emt? 2x sit Posty om, VIMUI2Ss9hnwy WMH .1veqs, 


“n't 26 umitosh sels Uf hy @tenitpoxehna ane . Lrtovee 


n@tcuiuqeq tihdins ube to aeeeis a) Shvsqed eaeeetpotg, 


talsne ,notteluqeg beitbabeadony -968T00R Iiiw some 


Sesumtsuasewan yitipiie Ide =e quote aps. €i~Ni eg ak 
P.O #3 4. vend seeidiing epiel alt ee Sesuqgmeg as) 


“cs ed* .25ey pulbesoous dose of ot0 tee |. eaemeelie, 


Ni Stemivorgin yiseelo opel) Itie quem ows Siefee2 $8 


O96819 doze baw “ode vontod Gas .«Qioup ane. Ants we notiasuqog o'z3. 
| y 
vi? mi 9¢bmy Seonebiy Se m7? 02 e¢iueay, ites pains wie2zid AL, 


ay -hetaaimels: of yleviesosnoxe Elin aod niu bt -o$ 


Js ’ 
7692 sSitnGioul: Setdas \, aiicnes Lares tyerx) ode Onininexs 


*ei0i Pidetors. Yieeeed et gads a2 edistive on al oxety . 


-bemhLanolegs ed biyata® 


Table 3 indicates that marriages are slightly 
overestimated. This can be explained by the fact 
that any individual of marital status "other" is, 

in the model, eligible for marriage. We recall 

that the state-variable marital-status can attain 
three state-codes, i.e. single, married, and other. 
The "other" includes the widowed, separated and 
divorced. Including the separated into the eligible 
population for marriage clearly ae enced positive 


biases. 


The divorce process results are given in table 
4, and are much better than expected. It is known 
that Canada passed through a transient period during 
the time of the simulation, insofar as the incidence 
of divorces is concerned, due to a change in the divorce 
law in 1969. This made the simulation of divorce quite 
difficult, and hence the results are on the whole quite 
pleasing. It is obvious, however, that the regionalization 
of the divorce probabilities should be seriously considered. 
In 1970, for example, the model significantly overestimates 
the number of divorces in Quebec while underestimating 
the number in Ontario. It is clear that the number of 
divorces in these two processes cannot be considered to 


be outcomes of the same stochastic, process. 


It can be seen from table 5 that the death process 
results in an underestimate in the number of simulated 
deaths. This underestimate is 9.4%, 13.6% 11.1% and 8.2% 
respectively in the years 1968 through 1971. Most of 
this error (approximately 8%) is a consequence of the 
initial population underestimate, while simulation 


error can account for another -4%. 
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Table 6 compares the sex-region populations produced by 
POLSIM with the same distributions as reported by Vital 
Statistics. It can be seen that there is a general 
underestimation on the part of the model. This can be 
explained by the fact that the Survey of Consumer 

Finance underestimates the total population (by excluding 
the Yukon and N.W.T., military personnel, and persons 

in institutions) and because the model itself underestimated 
the number of births in each of the simulated years. 

Over the four year period the underestimation in population 
that can be attributed to the model itself is 128,276 

(sum of underestimates in births less sum of underestimates 
in deaths). The total underestimate in the 1971 simulated 
population is seen from table 6 to be 780,260. Of this 
total, 16% can be crudely attributed to the error 

generated by the model while 84% can be attributed to 

the error in the initial year population. This is in 

fact an upper estimate of the model error. The fact 

that the initial population is too small to begin with 
implies that one would expect an underestimate in the 
number of births and deaths. Therefore part of the 
underestimate arising from the model is in fact attributable 


to the error in initial year population. 


The results of the labor force simulations are summarized 
briefly in table 7 and Figure 2. Table 7 presents the 
unemployment rates, by sex, that the model produces 


over the entire four year period and compares these 
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with the actual rates for the same period as measured 

by the labor force survey. Figure 2 plots the monthly 
aggregate simulated unemployment rate and the monthly 
aggregate actual unemployment rate over the four year 
period. It can be seen that the model tracks the 
unemployment rate very well, and there is no tendency 
HOrett «eo get “orf=track”™ as time procresere There is 

a slight tendency, however, for the simulated rate to 

be too high. This was expected, because the simulation 
parameters are adjusted to fit the higher unemployment 
Mates of the period April 1972 to April 1973. 9 This 
adjustment was such as to increase the resulting simulated 
unemployment rate slightly from that which would have 
resulted from the original equations. Since the original 
regression equations of the labor-force model had been 
fitted to data from the 1959-1969 period, the adjustment 
would be expected to simulate too many unemployed 

persons over the four years 1967-71. The adjusted 
equations would be expected to perform better over a 


four year period beginning in 1971. 


The Market Income simulation is summarized in 
Tables 8-10. For each of the component incomes (employment 
income, property income, and retirement income) distributions 
for the four simulated years are presented. As standards 
of comparison, the same distributions from the 1967 and 
1971 SCF surveys are also given. It can be seen that 


the simulations perform as one would expect; there is a 
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general tendency for the distributions to shift to the 
right as time progresses. Since no adequate income 
data exists for the years between 1967 and 1971, it is 
not possible to examine how well the simulations perform 
year by year. All that can be compared are the final 
results (the distributions for the year 1971). It 
should be noted that the comparisons for property 
gncome in Table 9 are not really meaningful, due to the 
fact that the 1967 survey did not distinguish between 
Gividend and interest income. As a result, the model 
simulated total property income with interest income 
transition matrices, and hence the resulting final 
simulated property income is not strictly comparable 


with the results obtained from the 1971 SCF survey. 
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Table 1 


Comparison of simulated and actual data* 
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Table 3 


Comparison of simulated and actual data* 
Marriages by Region: 1968-71 


Atlantic Quebec Ontario Prairies B.C. Canada 

1968 Simulated 18, 200 Sys (S10, 70,775 SB) SVAS) 20,850 199,400 
Actual 16,665 46,004 62,109 29,678 16,914 171,766 

1969 Simulated 18,500 SYS. 0/5) TALS DS) Sy i) 19,700 187,450 
Actual 17, 420 47,545 67,150 31, 378 18,284 182,183 

1970 Simulated 19, 800 55, 600 L302 5 34,900 19,950 ZOSSS/ 5 
Actual 17,875 49,607 68, 874 31,610 20,026 188,429 

1971 Simulated 19,950 56,475 74,350 35,000 21,475 207,250 
Actual 18,678 49,695 69,590 32,554 20,389 Ons 24: 

Table 4 


Comparison of simulated and actual data* 


Divorces by Region: 1968-71 


Atlantic Quebec Ontario Prairies BeGe Canada 

1968 Simulated ILA tey/S) Se 87/5 1025 3532) 750 20,550 
Actual 675 606 5,036 25105 2220 115343 

1969 Simulated 550 6,075 7, 800 3050 1,825 20,800 
Actual 1, 362 2,930 Ti o4S 5,648 MW iespa gle: 26,079 

1970 Simulated 1,400 6, 450 8,350 2,950 2,600 2s. 750 
Actual 1,414 4,865 125454 5,876 Sys UL ZO fap 

1971 Simulated 2, 000 5,975 8,450 2,800 2,400 21,022 
Actual 13423 53195 129 5,835 4,942 29,626 


* Source: Vital Statistics - Statistics Canada Catalogue 84-201. 
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Table 5 


Comparison of simulated and actual data* 
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Deaths by Region: 1968-71 
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Comparison of simulated and actual data* 
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Population Distribution by Sex and Region: 1968-71 
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19, 768, 300 
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43,400 
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Final ‘70 


20,244,100 
97.94 


53,800 
- 26 


56,100 
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44,400 
21 
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109, 700 


20, 669, 600 
100.00 


Final *71 


20,474,100 
97.88 


57,050 
27 


57,200 
~27 


46,300 
22 


48,100 
23 


112,000 
-54 


59,900 
«29 


27,800 
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£55700 
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7,800 
-04 


20,916,750 
100.00 
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73,450 
34 


42,200 
20 
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APPENDIX 


FLOW CHART OF THE POLSIM MODEL 
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THE IMMIGRATION BLOCK 


\ 
START \ 


READ ALL 


CONSTANT 
PARAMETERS 


TOTAL IMGRANTS. READ ALL 
UNEMPLOY. RATE —- =—_- =— VARIABLE 
WEIGHT FACTOR PARAMETERS 


DETERMINE THE 
TOTAL NUMBER OF 
PEOPLE IN EACH 
AGE-—SEX—PROV. 
MSTAT CATEGORY 


INSURE THAT NO. 
OF MARRIED MEN = 
NO. OF MARRIED 
WOMEN IN EACH 
PROVINCE 


INITIALIZE 
FAMILY UNIT 
NUMBER 


ENTER 
LOOP FOR 
ALL 10 
PROVINCES 


ENTER 
LOOP FOR ALL 
SINGLE MEN IN 

THE GIVEN 
PROVINCE 


CALCULATE 
THE STATE 
VECTOR FOR THE 
GIVEN SINGLE 
MAN 


YES 


Psu a 
L ‘ 


1°. Seared 
Ti LMA 


Wale 


oi .A 
. e22wivorre 
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gen, 


aa peta 
) 4 Bt WM) 
VivViAM ZUM 


: UavVieSVT 
SVL 


Sfa seule 


VECTOR OF ALL 
IMMIGRANTS IN 

PROVINCE ON 
OUTPUT TAPE 


WRITE STATE 


ENTER 
LOOP FOR 
ALL SINGLE 
WOMEN IN THE 
PROVINCE 


CALCULATE THE 
STATE VECTOR 
FOR THE GIVEN 
SINGLE WOMAN 


SINGLE 
WOMAN IN 
THE PROV 


MES 


ENTER 

LOOP FOR 
ALL FAMILIES 
IN THE PROV. 


CALCULATE 
THE HUSBAND'S 
STATE VECTOR 


CALCULATE 
THE WIFE’S 
STATE VECTOR 


FAMILY SET FAMILY 
HAVE SIZE FOR 
CHILDREN CHILDLESS 


? COUPLES 


DETERMINE 
HOW MANY 
CHILDREN 


SET FAMILY 
SIZE FOR WHOLE 
FAMILY 


CALCULATE 
THE STATE 
VECTOR OF ALL 
THE CHILDREN 


LAST 
FAMILY IN 
PROV. 
? 


NO 
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DEMOGRAPHIC BLOCK FLOW CHART 
(PHASE 1) 


—_ READ A FAMILY UNIT 
BASE YEAR 
FILE TARE 


EMIGRATION & DEATH 
REGISTRY MAGN, TAPE 


PICK 1ST IND/L 


NO RECORD !IND/L ON 
ve DEATH REGISTRY 


| 
| 
| 
we 


LAST 
PICK NEXT IND/L € IND/L 
? 


lYES 


HAS 
ANY BODY 
SURVIVED 
iY 


BES 


NO, 


PICK 1ST IND/L 


~ <> 
YES 


WILL 
SHE HAVE 
A BABY? 


CREATE A CHILD 


MARRIED 
MOTHER? 


NO 


MAKE HER AND BABY A 


NEW FAMILY UNIT 


} NO _ PICK NEXT IND/L 


YES 


(CONTINUED) 


jal va? 0a 


a0 | 
aVAR Tem 
\\eAn “ 


¢ 
a} 


Ova Teale) 


& Tue hia tae 


aT 


AG 


DECIDE 
WHO 
CONSIDERS 
DIVORCE. 


WILL 
THEY 
DIVORCE? 


CHANGE MARITAL STATUS. 


SPLIT FAMILY IN TWO. 
(CHILDREN WITH MOTHER) 


NO 


PICK 1ST NON-MARRIED IND/L. 


WILL 
IND/L 
MARRY? 


ulster MARK HIM FOR MARRIAGE. 


DECLARE HIM HEAD OF 
NEW FAMILY UNIT. 


PICK NEXT IND/L. 


YES 


PICK 1ST "NEW" FAMILY UNIT. 


E BASE FILE” 
AGN. TAPE 


WILL 
“NEW” FAMILY \yves| RECORD ‘'NEW” FAMILY ON 


EMIGRATE? EMIGRATION REGISTRY 


EMIGRATION & DEATH 
REGISTRY MAGN. TAPE 


NO 


“MARRIAGE POOL” 
MAGN. TAPE 


1S 
HEAD 

GETTING 

MARRIED 


RECORD FAMILY UNIT ON 
MARRIAGE POOL TAPE. 


RECORD FAMILY UNIT ON 
“AGED” BASE FILE. 


PICK NEXT NO 
“NEW” FAMILY UNIT. 


IYES 


FINISHED 
WITH BASE 
YEAR TAPE 


ese STOP 


DEMOGRAPHIC BLOCK FLOW CHART 
(PHASE 2) 


START 


ASSIGN Tin TO MARRIAGE POOL TAPE 


AND Ts TO A SCRATCH TAPE 


INITIALIZE PROVINCE r.1 
; READ A RECORD 


HAS 
DEPENDENTS 
? 


READ DEPENDENTS’ RECORDS 


(IGNORE HiM) 


? 


STORE DEPENDENTS IN MEMORY 
AND SAVE ADDRESS IN HEAD'S RECORD 


DETERMINE SEX, AGE, EDUCATION 
AND KEEP RECORD OF HIM/HER 


STORE HEAD’S RECORD ON DISK-FILE 
AND SAVE HIS/HER ADDRESS 
IN ORDER TO RECALL HIM 
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Using N(s,a,e) population of heads, by sex, age and 
education, prepare a list, L(s,i), which indicates that 
the i-th head of sex s will marry the L(s,i}-th 
head of opposite sex. If L(s,i) » O it means that 
there was not available mate for head (s,i). 


DETERMINE SEX s’ WITH MOST HEADS 
AND SEX s‘' WITH LEAST HEADS 


INITIALIZE i+ 1 
“HEAD (s‘,i) DID NOT FIND MATE” 
YES 
RECALL HEAD (s‘,i) FROM DISK a 


| NO 
JOIN HIM WITH DEPENDENTS (IF ANY) 
©) RECALL HEAD L(s',i) OF SEX s” 
| | RECORD THEM ON TAPE T, P & 


JOIN HIM WITH DEPENDENTS (IF ANY) 


© RECALL HEAD i OF SEX s’ 


JOIN HIM WITH DEPENDENTS (IF ANY) 


RESTRUCTURE HEADS AND DEPENDENTS TO 
FORM A FAMILY AND UPDATE THE MARITAL 
STATUS CODE. 
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UPDATE PROVINCE OF RESIDENCE 


FOR THE FAMILY 


RECORD FAMILY ON NEW BASE YEAR FILE Tout 


FINISHED ‘ 


YES 


IS 


YES 
REWIND Tin 
NO 
: STOP 


MES “i.e, consider the leftovers 
from all provinces.” 


ASSIGN Tin TO THE SCRATCH TAPE 
REWIND Tin 
a. =e ee se Biss: sa Chensee 


e 
: : “If there are any final leftovers they must 
ASSIGN T, TO BASE YEAR FILE Tout be simply recorded on the base year file.”’ 


“'Go back to process leftovers considered as 
if they are from dummy province + 11. 
Read from Tin - scratch tape and record 
results on Tout, the new base year file.” 
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ACTIVITY STATUS BLOCK 


READ INPUT 
PARAMETERS 


CALCULATE 
TRANSITION 
PROBABILITIES 


READ AN 
INDIVIDUAL 
RECORD 


ENTER LFCHG 


| THIS IS THE MAIN | 
SUBROUTINE. 


PERSON IS CHANNELED | 
DEPENDING ON EMPLOYMENT 
CATEGORY AND AGE. 


| MONTHLY SIMULATION IS | 
CARRIED OUT IF NECESSARY. 


STATE VARIABLES ARE 
| UPDATED. 


EXIT LFCHG 


NO LAST 
RECORD? 


MACRO FLOW CHART OF THE 
ACTIVITY STATUS BLOCK 
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THE ACTIVITY BLOCK 


ACTIVITY 
BLOCK 
DATA 
TAPE 


READ ALL 
CONSTANT 
PARAMETERS 


YEAR, NO. OF 


READ ALL 
VARIABLE MONTHS, INITIAL 
PARAMETERS Sa MONTH, 


UNEM. RATES 


CALCULATE 
LABOR FORCE 
TRANSITION 
PROBABILITIES 


DISAGGREGATE 
TRANSITION 
MATRICES TO 

MARITAL STATUS 

$ REGION 


READ AN 
INDIVIDUAL , 
RECORD 


YRACT — 23 
PLACE PERSON af AGE <14 WKSCHL — 40 
IN GRADE 9 ? WKNLF — 12 
THIS YEAR WKEMP — 0 


WKUNEM — 0 


| >14 


ACTIVITY IN 
FIRST MONTH 


(APRIL) IS 
“YRACT” 


TYPE= 
TYPE/10 
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PERSON 
RETIRES 
TYPE = 40 


PERSON WKEMP = 0 
RETIRED YES WKSCHL =0 

? WKUNEM = 9 ice 
WKNLF = 52 


YRACT = 22 


PERSON 

NO RETIRES 
NO 

PENSION 


PERSON 
RETIRES 
TYPE = 50 


WILL 
PERSON 
RETIRE? 


NEVER WKEMP =52 
UNEMPLOYED WKUNEM =0 
WKNLF =0 
WKSCHL =0 
YRACT = 20 
UNEMPLOYMENT 
POSSIBLE 


PERSON DIDN'T 
WORK LAST YEAR 


YES .. MUST ASSIGN 


INITIAL INCOME THIS 
YEAR TYPE = 
TYPE 10 


CONVERT LAST 
YEAR’S WAGES TO 
WEEKLY RATE 
WAGES = WAGES 
+ WKEMP 
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BLOCK DETERMINES 
TRANSITION FROM 
INPUT STATE IN INPUT 
MONTH (APR—&MAR) 
TO OUTPUT STATE IN 
OUTPUT MONTH 
(MAY—®APR) 


AUG 
ei 


NLF 


STATE 


MONTH AFTER 


INPUT STATE 
FOR TRANSITION 
PURPOSES IS 


TEMPORARY 
OUTPUT STATE 
= INPUT 


LOOP FOR 
INPUT MONTHS 


IN SCHOOL IN 
INPUT MONT 


DETERMINE 
INPUT MONTH 


FOR TRANSITION 


(1 = 120) 


FIND TEMP. 


OUTPUT STATE 


WILL BE EMP, 
UNEM, OR 
NLF 


OUTPUT 
MONTH = SEPT 
? 


STATE IN OUTPUT 


MONTH = 


TEMPORARY 
OUTPUT STATE 


TEMPORARY 
NO OUTPU~ STATE 
= INPUT 
STATE 


WES 


TEMPORARY 
OUTPUT STATE 
=NLF 


INPUT 
STATE = 
SCHOOL 
i 


YES 


YES 


APR OR MAY 


NO 


PERSON 
STARTING 
SCHOOL? 


YES 


OUTPUT STATE 


YES 


TEMPORARY DETERMINE 
OUTPUT STATE WHERE IN 
= NEW EDUCATION SCHOOL UPDATE 
EDUCATION 


STATE 
STATUS 


NO 
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DETERMINE 
WEEKS IN EMP, 
UNEM, NLF, 
AND SCHOOL 


YRACT = LAST 
MONTH'S (APRIL'S) 
OUTPUT STATE 


IN LABOR 
FORCE 
2? 


DOES 
MALE 
PERSON HAVE 
UNIVERSITY 

DUCATIO 


TYPE = 200 

(ALL FEMALES & 
ALL MALES WITH- 
OUT UNIVERSITY) 


TYPE = 100 
(ALL MALES WITH 
UNIVERSITY) 


PENSIONABLE 


EMPLOYMENT 
? 
NO 
TYPE = TYPE TYPE — re 
+40 +50 


ACTIVITY 


UPDATE 
TAPE 


WRITE 


OUTPUT aaa 
RECORD 


LAST 
INDIVIDUAL 
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THE MARKET INCOME BLOCK 


START 


MARKET 

INCOME 
PARAMETER 
TAPE 


READ ALL 
armre CONSTANT 
PARAMETERS 


KYR & 
RINFL 
ON CARDS 


READ ALL 
VARIABLE 
PARAMETERS 


CALCULATE 
CONSUMER 
PRICE INDEX 
VECTOR 


OUTPUT 
TAPE FROM 
ACTIVITY 
BLOCK 


READ AN 
ae INDIVIDUAL 
RECORD 


RECORD OUTPUT 
RECORD 


YES 


WRITE OUTPUT 
RECORD & EXIT 


140, 150 


y 24,25 40 ; 50 
240, 250 3 14, 15 4.5 2 


INCREMENT “OTHER” 
MONEY INCOME 
BY ACPI 
SUM TO GET 


TOTAL INCOME 


TYPE DEFINITION 


OUTPUT 
TAPE FROM 
MARKET 


140, 150 EMPLOYED: FIRST YEAR IN 
240, 250 LABOR FORCE 


WRITE THE 
OUTPUT STATE _ 
VECTOR 


STUDENT OR NON LABOR FORCE 


EMPLOYED CLASS A PERSON 


RETIRED 
LAST NO 


24, 25 EMPLOYED CLASS 8 PERSON RECORD 


40 RETIRED FIRST YEAR NO PENSION 


YES 
RETIRED FIRST YEAR 
ELIGIBLE FOR PENSION 
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INITIAL AND STUDENT INCOMES 


RETIREMENT 
AND EMPLOYMENT 
INCOME BOTH 
=0 


RETIREMENT 
INCOME = 0 


FIND MEAN 1971 
FULL EMPLOYMENT 
INCOME (MFE}!) 
IN 1971 DOLLARS 


FIND KYR 
MFEI IN 
1971 DOLLARS: 
REAL GROWTH 


MONEY GROWTH: 
FIND KYR MFEI 

IN KYR DOLLARS 
=MFEI (KYR) 


TYPE 
OF 
PERSON 
? 


EMPINC = 


STUDENT 


EMPINC = 
MFEI (KYR) 


TYPE = TYPE/10 


CLASS B 


(DEPENDS ON 
AGE, SEX, EDUCATION 
MSTAT, PROVINCE) 


(GROWTH FACTOR 
DEPENDS ON AGE, 
SEX, EDUCATION) 


(USES CHANGE 
IN CPI) 


EMPINC = 


52 


MEENKYR) , WKEMP 


- 
= Jv inlin | 


(Ral fae 
. | nt 
; | 


gaser = 307 


(DEPENDS DETERMINE 
ON SEX AND INITIAL 
LAST YEARS RETIREMENT 

INCOME) INCOME 


EMPLOYMENT 
INCOME =0 
TYPE = TYPE/10 


RETIREMENT 
INCOME = 


CONVERT PRESENT 
ANNUAL WAGE 
TO 1970 DOLLARS. 
FACTOR = A CPI 


FIND PRESENT 
INCOME CLASS AND 
RELATIVE POSITION 

IN THAT CLASS 


DETERMINE 
ANNUAL INCOME 
TRANSITION. NEW 
INCOME CLASS !S 
IN 1970 DOLLARS 


INCOME IN NEW 
INCOME CLASS IS 
IN SAME RELATIVE 
POSITION AS IN 
OLD INCOME CLASS 


CONVERT NEW 
INCOME TO 
KYR=1 DOLLARS 
FACTOR = ACPI 


EMPLOYMENT 
INCOME =0 
RETIREMENT 
INCOME = 0 
TYPE = TYPE/10 


RETIREMENT 
INCOME =0 


CONVERT PRESENT 
WEEKLY WAGE 
TO 1969 DOLLARS. 
FACTOR = ACPI 


FIND PRESENT 
WAGE CLASS AND 
RELATIVE POSITION 
IN THAT CLASS 


DETERMINE 
(TRANSITIONS WEEKLY WAGE 
DEPEND ON TRANSITION. NEW 
AGE AND SEX) WAGE CLASS IS 
IN 1969 DOLLARS 


WAGE iN NEW 
WAGE CLASS IS 
1 iN SAME RELATIVE 
POSITION AS IN 
OLD WAGE CLASS 


CONVERT NEW 
WAGE TO 
KRY=1 DOLLARS 
FACTOR = A CPI 


EMPLOYMENT 
INCOME = NEW 

} WAGE RATE X 

} WEEKS EMPLOYED 
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PROPERTY INCOMES 


DETERMINE 
AGE CLASS AND 
TOTAL INCOME 
CLASS 


RANDOMLY 
DETERMINE 
NEW DIVIDENDS 
TRANSITION MATRIX 
DEPENDS ON AGE 
AND INCOME 


INITIAL 
DIVIDENDS > 
$250 
? 


RANDOMLY 
DETERMINE 
NEW DIVIDENDS 
TRANSITION MATRIX 
DEPENDS ON AGE 
ONLY 


RANDOMLY 
DETERMINE 
NEW INTEREST 
TRANSITION MATRIX 
DEPENDS ON AGE 
AND INCOME 


INITIAL 
INTEREST 2 
$250 
rd 


RANDOMLY 
DETERMINE 
NEW INTEREST 
TRANSITION MATRIX 
DEPENDS ON AGE 
ONLY 


